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ABSTRACT 

We discuss the potential of using the He I 584 A forest to detect and study He n reionization. 
Significant 584 A absorption is expected from intergalactic He n regions, whereas there should be 
no detectable absorption from low density gas in He ill regions. Unlike He II Lya absorption (the 
subject of much recent study), the difficulty with using this transition to study He n reionization 
is not saturation but rather that the absorption is weak. The Gunn-Peterson optical depth for this 
transition is r ~ O.lxHeii A 2 [(1 + z)/5] 9 / 2 , where ^hoII is the fraction of helium in He n and A is 
the density in units of the cosmic mean. In addition, He I 584 A absorption is contaminated by lower 
redshift H I Lya absorption with a comparable flux decrement. We estimate the requirements for a 
definitive detection of redshifted He I absorption from low density gas (A « 1), which would indicate 
that He II reionization was occurring. We find that this objective can be accomplished (using coeval 
H I Lya absorption to mask dense regions and in cross correlation) with a spectral resolution of 10 
and a signal-to-noise ratio per resolution element of ~ 10. Such specifications may be achievable on 
a few known z ~ 3.5 quasar sightlines with the Cosmic Origins Spectrograph on the Hubble Space 
Telescope. We also discuss how He I absorption can be used to measure the hardness of the ionizing 
background above 13.6 eV. 

Subject headings: cosmology: theory - cosmology: intergalactic medium - quasars: absorption lines 



1. INTRODUCTION 

Starlight produced by the first galaxies is the lead- 
ing candidate for ionizing the hydrogen as well as singly 
ionizing the helium at z ~ 6. It takes a harder 
source of radiation to doubly ionize the helium, so 
the reionization of this species is likely deferred un- 
til z ~ 3 w hen quasars p roduce a sufficient hard UV 
background (iMadau et al] 11999b IFurlanetto fc Ohl 120081; 
IMcQuinn et all 12009). However, the helium could have 
been doubly ionized at nearly the same cosmic time 
that hydrogen was reionized if more exotic sources ion- 
ized the hy drogen, such as the first generation of metal- 
free stars (iBromm et al.l 12001b IVenkatesan et al.l 120031: 
Tumlinson et al. _ 2004f) or min iquasars ( Madau et al.l 



ization (IFurlanetto fc Ohl [20(51 IMcQuinn et all [2009; 
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1999b iVolonteri fc Gned in 2009). In a third potential 



scenario, early sources doubly ionized the helium and 
then shut off. Afterward, the He II reco mbined such that 
quasars could again reionize it at z ~ 3 (|Wvithe fc Loebl 
2003i; IVenkatesan et al.ll2003l ). 

If He II reionization were completing at z ~ 3, an 
epoch for which there are numerous observations of 
the intergalactic medium (IGM), it should be an eas- 
ier task to definitively detect this process compared 
to detecting z > 6 reionization processes. Further- 
more, if He II reionization were ending at z ~ 3, it 
should have significantly affected the temperature of 
the intergalactic gas and the ultraviolet radiation back- 
ground. These motivations, along with recent addi- 
tions to the Hubble Space Telescope (HST), have in- 
spired a significant effort of late to understand the sig- 
natures and the detection prospects of He n reion- 
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Three separate observations of the z ~ 3 IGM suggest 
that He n reionization was ending around this redshift: 
First, several studies have measured the temperature of 
the intergalactic gas from the widths of the narrowest 
lines in the H I Lya forest, and the majority of these 
studies have found evidence for an increase in the IGM 
temperature of ~ 10 4 K between z ~ 4 and z ~ 3 , 
before a decline to lower redshift (|Schave et al.l 120001 : 
iRicotti et all 120001: ILidz et al.ll2009ft . These trends have 
been attributed to the heating from He II reionization. 
Second, observations of He II Lya absorption from gas 
at 2.8 < z < 3.3 show tens of como ving Mpc (cMpc) re- 
gions with no det ected transmission (jReimers et al.lll997t 
iHeap et al.ll2000l ). which may indic ate that He II r eion- 
ization was not comple te. Thirdly, iSongailal (119981 ) and 
lAgafonova et al.1 (|2007f ) detected evolution in the column 
density ratios of certain highly ionized metals at z rs 3, 
which they argued was due to a hardening in the ioniz- 
ing background around 50 eV and, thus, the end of He n 
reionization. 

However, the interpretations of all of these indications 
for He II reionization are controversial. Temperature 
measurements of the IGM are difficult, and not all mea- 
surements detected the aforementioned trends. It is of- 
ten argued that He II Lya absorption saturates at He II 
fractions that are too small (~ 10~ 3 at the mean den- 
sity) to study He II reionization (although, see McQuinr] 
2009). Lastly, inferences from metal lines require signif- 
icant modeling, and studies have reached different con- 
clusions regarding the d egree of their evolution at z ~ 3 
(Bok senberg et all 12001) . 

This paper discusses intergalactic absorption by the 
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Is 2 : S — lsnp 1 P° transitions of He I as an unsaturated 
observable of He n reionization. We primarily focus on 
the strongest and longest wavelength of these absorption 
lines, the He I Is 2 X S - ls2p 1 P°, 584 A line. For a 
given optical depth in the H I Lya forest, the amount of 
absorption in the He I forest can be directly estimated 
if both the fraction of helium that is He n and the ratio 
of the H I and He I photoionization rates are known. In 
addition, because the He I ionization edge is relatively 
close to that of hydrogen, the photoionization rate that 
He I experiences is similar to this for hydrogen. Thus, 
the most important determinant of the amount of He I 
584 A absorption is the He n fraction. 

Other studies have discussed the He I forest, but did 
not focu s on its usefulness as a probe of helium reion- 
ization. iTripp et al.l (1990) was the first to discuss in- 
tergalactic absorption by the He I 584 A, and this study 
attempted unsuccessfully to det ect this absorption in the 
spectrum of a z — 1.7 quasar. iReimers fc Vogell (1993) 
targeted He I absorption from z ss 2 H I Lyman-limit 
systems, and reported the first (and at present only) de- 
tection of intergalactic He I 584 A absorption. Finally, 
iMiralda-Escude fc Ostrikerl (|1992h showed that the He I 
forest could be a useful probe of the hardness of the ul- 
traviolet background, and they focused on using this ab- 
sorption to rule o ut a particular n onsta ndard model for 
the dark matter. iSantos fc Loebl (|2003f ) followed up on 
this idea, arguing that the z ~ 5 He I forest could be a 
useful diagnostic of the hardness of the ultraviolet back- 
ground after hydrogen reionization. 

As with the He II Lya transition at 304 A, the He I 
584 A transition falls blueward of the hydrogen Lyman 
limit and, therefore, is subject to continuum absorp- 
tion by hydrogen, in this case from H I systems with 
z > z' = 3.2(1 + ZHei)/5 — 1. This continuum absorp- 
tion may even be worse in terms of obscuring the He I 
584 A forest compared to the He II Lya forest because 
this spectral region is more affected by higher redshift H I 
continuum absorbers. It is unlikely that there are more 
than a handful of quasar sightlines at z > 4 with suffi- 
cient near ultraviolet (NUV) flux for detection with the 
present generation of instruments (and almost certainly 
not a t z > 6, during He I reionization; Santos & Loeb 
2003). Therefore, our focus is on applications of He I 
absorption at z < 4. 

We have performed a cursory search for candidate He I 
sightlines among the relatively few published Hell Lya 
forest spectra that extend redward to 584 (1 + zqso) A. 
Of note, QSO OQ 172 (z = 3.54) has F x ~ 2 x 
10~ 16 erg s" 1 cm" 2 A^ 1 at 584 (1 + z) A (jLvons et al.l 
\EM), and QSO 0055-269 (z = 3.67) has F A «lx 10" 16 
(Gabor Worseck, priv. com.). In fact, we were surprised 
to find that most He n sightlines in our search had sig- 
nificant flux in the relevant band. 

Many of the existing He II sightlines had been selected 
by their far ultraviolet flux. NUV selection would be 
more optimal for identifying candidate He I forest sight- 
lines. For example, HS 1140 +3508 (z = 3.15) is ob- 
scured in the far ultraviolet by a Lyman-limit system, 
but has a NUV flux of F x ~ 2 x 10~ 16 (Gabor Worseck, 
priv. com.). 

Another difficulty with the He I forest is that line ab- 



sorption from foreground systems can contaminate the 
He I forest, the most important of which is H I Lya 
absorption. H I Lya absorption from a system with a 
redshift of zhi,1216 = 2.4(1 + ZHei,584)/5 — 1 falls di- 
rectly on top of the He I absorption from a redshift of 
ZHei,584- Fortunately, the Lya forest is quite thin at rel- 
evant zhi,i216 (with a flux decrement of several percent), 
and we find that both this contaminant's mean transmis- 
sion and variance tend to be comparable to that of the 
He I forest. Also, the He I forest correlates strongly with 
the coeval H I Lya forest, which allows it to be extracted 
in cross correlation despite this contamination. 

This study is timely because the HST reservicing mis- 
sion installed the Cosmic Origins Spectrograph (COS) 
in May 2009. COS is capable of measurements of Hcl 
584 A absorption at redshifts relevant to He II reioniza- 
tion (2.8 < z < 4.5), and ground-based spectrographs 
can cover higher redshifts {z > 4.3, corresponding to 
> 3100 A). COS is able to achieve higher signal-to- 
noise ratios than previous instruments in the ultravio- 
let. A 60 hr exposure with COS would achieve a signal- 
to-noise ratio of 10 at R ~ 30, 000 for a flux of F\ = 
2 x 10~ 16 erg s" 1 cm" 2 A" 1 (the flux of QSO OQ 172 
and HS 1140 +3508) at 2600 A (He I absorption from 
z = 3.5). A 40 hr exposure with COS would achieve a 
signal-to-noise ratio of 10 for this flux at R « 20, 000 at 
2300 A (z = 2.9)0 

This paper is organized as follows. Section [5] discusses 
the physics of the He I forest. Section [3] contrasts sim- 
ulated spectra for this forest under different assump- 
tions regarding the ionization state of the helium. Sec- 
tion|4]quantifies the spectral quality an observation must 
achieve to verify whether He II reionization was occurring 
from the He I forest. Appendix IA1 outlines how to mea- 
sure the hardness of the ionizing background with He I 
absorption, and Appendix [B] derives formulae for the sig- 
nificance with which He I absorption can be detected in 
cross correlation with the coeval Lya forest. 

This paper assumes a flat ACDM cosmology with 
h = 0.7, n b = 0.046, n m = 0.28, a s = 0.82, n s = 1, 
and Yne = 0-24, con sistent with recent measurements 
(|Komatsu et al.l [20090 . However, the simulation used to 
calculat e the Lya forest spectra at z < 1.5, the D5 simu- 
lation in lSpringel fc Hern quist (2003), assumes a slightly 
different cosmology, with the most notable differences 
being ttf, — 0.04 and as = 0.9. The photoionization 
and recombination r ates used in this study are from 
IHui fc Gnedinl (fl99l . 

2. THE He I FOREST 

In photoionization equilibrium, the He I fraction is de- 
termined by the relation 



I Hcl 

Photoionization equilibrium is a good approximation be- 
cause the timescale to rea ch equilibrium (r^ T ~ T^ 1 
3 x 10 4 yr at z ~ 3; e.g., iFaucher-Giguere et al.lfeOOSaf ) 

3 http://etc.stsci.edu/webetc/. As of mid-January 2010, the 
dark current for the COS NUV gratings is 5 times higher than 
specifications (and than what is assumed by this ETC), and, thus, 
a significantly longer observation is required unless this issue is 
resolved. 
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was much shorter than all other relevant timescales. 
Here, a x , n x , x x , and T x = J^L a (dE/E) a x {E) J{E) 

are respectively the Case A recombination coefficient 
the number density, the ionization fraction, and the pho- 
toionization rate for ion X (or subscript e for electrons) , 
and E^ 11 and a x are the ionization potential and the 
photoionization cross section. Lastly, J(E) is the inci- 
dent specific intensity integrated over solid angle. We 
will sometimes imprecisely write XHeii = 1 for a gas par- 
cel even though a small fraction (~ 10 -5 A&) of the he- 
lium is in xhoI- 

If the density and Fuel are known (since n e is ef- 
fectively known in the IGM after hydrogen reioniza- 
tion up to density), a measurement of XHei can be 
used to determine a;Heii (equation [TJ . Fortunately, co- 
eval H I Lyman-series absorption provides an estimate 
for the density of a He I absorber. In addition, the 
He I photoionization rate (like this for the H i) is ex- 
pected to have been essentially spatially independent ow- 
ing to the long mean free path of He i-ionizing photons 
(luei) and the large number of sources in a volume of 
~ ijftj. It is estimated that luei > lm ~ 300 cMpc at 
z ~ 3 (iMadau et all 1 19991: IFaucher-Giguere et~aTll2008at 
Prochask a et al.l I2009T ) . Therefore, Tnei was effectively 
just a single number in all of the IGM. 

This paper focuses primarily on the longest wavelength 
allowed ground-state transition for He I, the Is 2 1 S — 
ls2p 1 P°, 584 A transition. Its oscillator strength is 3.8 
times larger than the next strongest transition of He I, 
the Is 2 X S - ls3p 1 P°, 537 A transition. He I 584 A 
absorption is observable in a suitable quasar absorption 
spectrum between the wavelengths of 584 (1 + zqso) A 
and 304 (1 + zqso) A, where the He II Lya forest begins. 
This wavelength range corresponds to He I absorption 
from 0.5 z QS o < z < zqso0 

To redshift across the 584 A transition, a photon ex- 
perienced the optical depth from gas in the Hubble flow 
of 
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(what is termed the He I Gunn-Peterson optical depth; 
iGunn fc Petersonlll965t iTripp et al.lll990l) where 
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is the H I Lya Gunn-Peterson optical depth, is the 
gas density in units of the cosmic mean, T4 is the tem- 
perature in units of 10 4 K, and the —0.7 exponent owes 
to the temperature dependence of the recombination co- 
efficient (cf. equation [1]) . 

4 Case A is most relevant for H I and He I because this paper's 
focus is on intergalactic systems in which these species are highly 
ionized. 

a He I continuum absorption by intervening systems will produce 
breaks in the spectrum starting blueward of 504 (1 + 2qso) A (the 
redshifted He I limit), but these discontinuities should be minimal 
(< 10%) for sightlines selected against having strong H I continuum 
absorption. 



Overdensities of a few and greater at z ~ 3 had de- 
coupled from the Hubble flow and were collapsing or had 
collapsed. In these regions, the Gunn-Peterson optical 
depth no longer describes the absorption. Instead, such 
regions appear as distinct absorption lines with widths 
of 10s of km s -1 , and their optical depth in He 1 584 A 
is 
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where Tjjt™ 2 i6 is the H 1 Lya optical depth, and we have 
approximated the line profile as a tophat with velocity 
width Av x (and 1 < A-uhi/AvhcI < 2, with 2 being the 
limit of pure thermal broadening). Systems with TVhi > 
10 14 cm -2 are common in the Lya forest. Each sightlinc 
inter sects — 100 such system s between 2 = 3 and z = 4 
(e.g- lPress fc Rvbicki|[l993l) . 

2.1. He I Photoionization Rate 

The He 1 584 A optical depth depends on the value 
of ThcI in addition to the field of interest, XHeii- The 
calculations in this paper assume ThoI = Thi, and 
Phi is chosen to match measurements of the H 1 Lya 
mean transmission (requiring r H i ~ 10~ 12 s _1 ; e.g., 
IFaucher-Giguere et al.l l2008af K This choice is moti- 
vated by ultraviolet backgrou nd models, which find 
r He i » T H i at z = 2 - 4 .5 (|Haardt fc Madaul 119961 : 
IFaucher-Giguere et al.| [2009) . The calculation of Thi and 
Tnei in ultraviolet background models requires two in- 
puts: (1) a model for the source emissivity at frequen- 
cies near the ionization potential for H 1 and He 1, and 
(2) the measured H 1 column density distribution, which 
determ ines the attenuation by H 1 cont inuum absorption. 

The IFaucher-Giguere et al.l (|2009f) model predicts 
ThcI = 0.8 — 1.0 Thi at 2 < z < 5. This model assumes 
that both stars and quasars contributed to the meta- 
galactic ultraviolet output, with an increasingly impor- 
tant stellar contribution with redshift at z > 2.5. This 
model uses the result of the stellar population synthesis 
code PEGASE for the spectral index of starlight between 
the He 1 and He 11 ionization edges. This code yields a 
similar value for stars (a ~ —1, defined as the pow er- 
law slope of the specific intensity; IKewlev et al.ll2001l) to 
what is expected for quasars (a ~ —1.6: iTelfer et al.l 
2002). These si milar a result in the near co nstancy 
of Tft e i(z) in the IFaucher-Giguere et al.l (|2009l ) model. 
Given these input a, H 1 continuum absorption then 
resulted in the spectral index of the average radiation 
background being hardened above 13.6 eV by roughly 
3 (/3 — 1) over the input a for < 2, where /3 is the 
power-law index of the H I column-density distribution . 
The background model of IFaucher-Giguere et al.l (2009) 
assumes /? « 1.5, resulting in a hardening of Aa « 1.5. 

However, there is modeling uncertainty in the input 
a for stars that ste ms from uncertainty in the spectrum 
of Wolf-Rayet stars (jKewlev et al.|[200ll ). in the average 
stellar metallicity, and in the stellar initial mass function. 
If the spectral index of stars were significantly softer, 
the ratio rHei/rm would be smaller than in the fiducial 
model in IFaucher-Giguere et al.l (|2009f) . In the extreme 
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example that the stellar contribution dominated the ion- 
izing background and had a spectral index of —3 instead 
of —1, ThoI would have been reduced relative to Thi by 
a factor of w 3. This boosts TH e i,584 by a factor of 3. 
In addition, there is u ncertainty in (3 at the ±0.3-level 
(Prochask a et al.ll2009f ) that affects Thci/Thi at the fac- 
tor of two-level. 

Previous calculations of the ionizing background in- 
cluded only H I continuum absorption to calculate (given 
a source emissivity model) the ionizing background be- 
tween and 0.75 E$£ u . These models ignored He I 
continuum absorption. This approximation may have 
affected their estimate for Tnei because the He I contin- 
uum absorption above E^^ is comparable to that of H I 
for optically thin systems. Namely, for gas of primordial 
composition in photoionization equilibrium 

« 0-5 zhcii ^ r^r , (5) 

t Hoi 

where r™ nt (E) is the continuum absorption of ion X at 
energy E. 

However, most of the absorption of background pho- 
tons with energy ss -EhS at z ~ 3 derived from systems 
with N m « 10 18 cm" 2 (for which rgft^f^) 1). 
These systems were much more self-shielded to H I 
Lyman-limit photons than to those at the the He I ion- 
ization edge. Thus, the effective value of Thi/ThcI in 
equation ([5]) was much lower in these dense absorbers, 
resulting in H I continuum absorption dominating over 
that of He I at k E££j. We estimate in calculations not 
included that the break at ES£i owing to He I continuum 
absorption in the spatially-averaged UV background at 
z ~ 3 amounted to only a couple percent. This justifies 
our use of a background model that does not include He I 
continuum absorption. 

3. CALCULATIONS 

Figures Q] and [2] show calculations of the He I forest 
at z — 4 and z = 3, respectively. These calculations 
use a randomly selected skewer of length rs 60 cMpc 
(Az « 0.1) through the 25 cMpc/fe, 2 x 512 3 particle 
SPH simulati on described inlLidz et al.l (120091) . The D5 
simulation in iSpringel fc Hernquistl (|2003| ) is also used 
for calculating the low-z Lya contamination^ The mean 
transmission in the H I Lya forest has been rescaled in 
these calculations by adjusting Thi to match the ob- 
served values. In par t icular, we use the Tni-values of 
iFaucher-Giguere et al.l ()2008bl ) for z > 2 and those in 
iKirkman etldl (|2007t) for z < 2. The He I 584 A for- 
est is calculated in post-processing with the conservative 
assumption that ThoI = Thi fSection l2.1j) . 

The three large panels in FiguresQ]and[2]makc different 
assumptions regarding the ionization state of the helium. 
The leftmost assumes that the He II fraction is given 
by XHeii = 1 everywhere, the middle that the helium is 
mostly He ill and that the fraction in He n versus He in 
is determined by a weak He n-ionizing background with 

6 The two simulations were used because the Lidz et al. (2009) 
simulation was not run to low enough redshifts to use to calculate 
the low-redshift Lya forest contamination. Both simulations were 
run to z = 2. We have compared calculations of the z ~ 3 absorp- 
tion in these two simulations (which have a mass resolution that 
differs by a factor of 10) to verify convergence. 



ThcH = 5 x 10~ 16 s _1 , and the rightmost assumes an even 
weaker background with ThcH = 10~ 16 s — 1 . The value 
ThcH = 5 x 10~ 16 s _1 is approximately the upper bound 
that has been set from observations of the most opaqu e 
regions in the He n Lya forest at z » 3 (|McQuindl20lM ). 

The densest absorbers may become self-shielded to 
He n-ionizing photons and can experience an even 
smaller ThcH- Self-shielding is included in these calcu- 
lations with a simplified 1-D algorithm for the radia- 
tive transfer. This algorithm assumes that regions with 
Ab < 1 experience the quoted THeii- For the denser re- 
gions, this algorithm locates segments along a given 1-D 
skewer that are bounded by At = 1. Each segment is 
composed of tens to hundreds of grid cells, each with 
width Ri 1 km s" 1 . 

For each overdense segment, our scheme attenuates the 
ionizing background (half of which is assumed to enter 
from each side) based on the amount and distribution 
of He II. The scheme starts with the assumption that 
ZHeii,i = 1 and Ji{E) = J(E) void exp[-a He ii(£) JV He n,i] 
for all cells across the segment, where i labels the cell 
number, and -/VH e ii,i and Ji{E) are respectively the He n 
column density and incident intensity to cell i from < 
i. Similarly, there is a contribution to Ji(E) from > 
i. Next, the scheme iterates to converge to rHcii,i and 
£Heii,i, assuming photoionization equilibrium. For our 
calculations, J(E) vold is set to have a spectral index of 
— 1.5, as expected for quasars. 

He II Lya forest sightlines show that THeii fluctu- 
ates wildly on scales of > 10 cMpc (|Zheng et al.l l2004t 
iFechner fe Reimersl 120071 ). We do not attempt to model 
these fluctuations here, but will comment on how they 
could affect our conclusions. These fluctuations will spa- 
tially modulate the number of dense self-shielding regions 
compared to the uniform-FHoii case. However, because 
the He I forest is sensitive to ZHeii ~ 1, ThcH fluctuations 
are less of a concern for studying He II reionization with 
the He I forest compared to with the He n Lya forest 
(which is sensitive to XHeii < 1). 

The top subpanels in Figures [T] and [2] show the value 
of XHeii that results from this algorithm, where Hubble's 
law is used to relate position to velocity. The blue curve 
is 4 x XHeii and the black curve is A;,. The value of 
iTHcii can be quite large for the weak backgrounds that 
are assumed. For A& = 1, T = 10 4 K, and ThcII = 
5 x 10" 16 s- 1 (1 x 10~ 16 s- 1 ), xhcII equals 0.06 (0.25) at 
z = 3. This number becomes 0.11 (0.39) at z = 4. 

The bottom subpanels in each larger panel in Fig- 
ures Q] and [2] zoom in on the residual transmission 
in the He II Lya forest (absorption at wavelengths of 
304 [1 + z] A). Note that in all of the cases there is min- 
imal transmission in this forest because the He II Lya 
transition saturates for xudi ~ 10 -3 A^ 1 . There is only 
a detectable amount of transmission in the case with 
r He ii = 5 x 10~ 16 s" 1 and z = 3 (middle panel, Fig.©. 
Because He n reionization was patchy, regions of trans- 
mission in the He n Lya forest can occur even if reioniza- 
tion was not complete. In contrast, He I absorption has 
the potential to reveal whether those opaque neighboring 
regions were in fact He n regions because it is sensitive 
tO XHoII ~ 1. 

The middle subpanels in Figures [1] and [2] show the 
transmission in the H I Lya forest (blue curves) and He I 
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Fig. 1. — Plot of different quantities from a simulated skewer through the IGM at z = 4 with length £3 60 cMpc (Az = 0.08). The three 
large panels make different assumptions about the ionization state of the He II. The leftmost panel assumes iCHell = 1; the middle that the 
He II is kept ionized by the background ThcH = 5 x 10 — 1B s — , and the right the same but ThoII = 10 -16 s _1 . For the latter two cases, 
self-shielding by dense regions is included with an approximate radiative transfer method that is described in the text. In each larger panel, 
the top subpanel shows Aj (black curve) and 4 X xjjell (blue curve), the second subpanel down shows the transmission in the He I 584 A 
(black curve) and H I 1216 A (blue curve) forests, and the third subpanel down zooms in on the He I 584 A transmission (black curve). 
The third panel also shows the observed transmission (blue curve), which, in addition to He I 584 A absorption, includes the contamination 
from foreground H I Lya absorption. The highlighted regions represent denser gas in which thi,i216 > 3 in the coeval Lyo forest. Lastly, 
the bottom subpanel zooms in on the transmission in the He II Lya forest. 
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Fig. 2. — Same as Figure [T] but at 



584 A forest (black curves). These panels demonstrate 
that the transmission in the H I and He I forests is highly 
correlated. The He I absorbers that correspond to the 
deepest H I Lya forest lines are the most visible, and the 
weaker He I lines (systems that are not dense enough to 
self-shield) disappear in the cases in which the helium is 
mostly doubly ionized. 

The third subpanel down in each larger panel in Fig- 
ures [T] and [2] zooms in on the the He I 584 A trans- 
mission (solid black curves). These subpanels also in- 
clude mock realizations of the observed spectrum at these 
wavelengths (i.e., He I 584 A plus foreground H I Lya ab- 



sorption; solid blue). Each panel uses a different skewer 
through the IGM to calculate this foreground absorption. 
The highlighted regions in these panels represent the lo- 
cations with thi,i216 > 3 for absorption coeval to that of 
the He I. Outside of the highlighted regions, the amount 
of He I absorption is substantially different between the 
three cases in each figure. Trace amounts of absorption 
remain in essentially just the xhoII = 1 case. A detec- 
tion of He I absorption in these regions would indicate 
that He n reionization was occurring. The next section 
quantifies the prospects for detecting this absorption. 
The top panel in Figure [3] plots the effective optical 
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depth (defined as r e g = — log T, where T is the aver- 
age transmission). This quantity is plotted for both the 
He I 584 A forest and the foreground Lya forest. The 
He i 584 A value of T e g is shown for the model with 
^Hcii = 1 (thick dashed black curve) and the model with 
ThoII = 10~ 16 s^ 1 (thick solid red curve). The value of 
r ff is comparable in these two models and comparable 
to that of the foreground Lya absorption (thick green 
dashed curve). Yet, r e g is 3 times larger at z = 3 in 
the XHcii = 1 model compared to the TroII = 10~ 16 s _1 
model, which could potentially allow these models to be 
distinguished. Also, the value of r e ff in the case with 
ThoII = 5 x 10~ 16 s _1 (not shown in Fig. [3]) is a factor 
of 10 smaller at z = 3 than the XHeii = 1 case. If THeii 
fluctuates spatially such with ThcII > 5 x 10~ 16 s _1 (ap- 
proximately the mi nimum value der ived from He II Lya 
analyses at z 3; McQuinn 2009), the signal will be 
even smaller than in this case. 

The thin curves that represent He I 584 A absorption 
in the top panel in Figure [3] are the same as the thick, 
except regions are masked with thi,i216 > 3 in the co- 
eval Lya forest (and the thin curves for the model with 
LhoII = 10~ 16 s" 1 are multiplied by 10). This figure also 
demonstrates that the transmission is much different be- 
tween the ^Heii = 1 and ThoII = 10~ 16 s _1 models if 
denser regions that have thi,i216 > 3 are masked. 

3.1. Recombining Regions and Double Reionization 
Scenarios 

We have shown that the amount of He I absorption is 
substantially different between the case where the He II 
was ionized by a weak background and the case xhcII = 1- 
There is also the exciting possibility that the He II was 
recombi ning after being ionized by an early generation of 
sources ( Venkatesan ct al. 2003) or after a nearby quasar 
had turned off. In this case, the He II fraction of a 
gas parcel is xhcII = 1 — cxp(— At/t roc ), where At is 
the time since the He n-ionizing background turned off, 
i rcc = (a]| cII n e ) _1 , and a]| cII is the Case B recombi- 
nation coefficient. At z — 3 and T = 20, 000 K, t rcc 
is roughly equal to 0.8 of the Hubble time at the mean 
density. 

Thus, the helium in underdense regions at z ~ 3 would 
have remained doubly ionized for At > H(z)^ 1 after 
the background turned off. The density dependence of 
a^Hoii is different in this case than in the other cases 
we have considered, such that there is the possibility 
that it also can be distinguished using He I absorption. 
However, it takes an extremely weak background with 
LhoII > a Hcii n e ~ 1.2 x 10~ 17 Ah s" 1 to coun teract re- 
combinations in a region. iMcQuinn et al.l ()2009l ) found in 
simulations that such a background developed soon after 
He II reionization by quasars was underway. Therefore, 
we consider the recombining scenario to be less likely 
than the others at 3 < z < 4, but nevertheless a tanta- 
lizing possibility. 

4. VERIFYING THAT He o REIONIZATION WAS 
OCCURRING 

The z ~ 3 He n reionization process is expected to 
have been very patchy. If quasars are the source of the 
ionization, models predict large-scale patches of primar- 
ily He II or He in gas that spanned many tens of cMpc 



foreground HI Lya 
same but masking x > 0.2 
Hel 584A, x Hen = 1 
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Fig. 3. — Effective optical depth (r ef f) and standard deviation 
in the transmission (ct) as a function of the redshift of the 584 A 
absorption (2584). These calculations assume Fuel = Thi, and Thi 
is chosen to match observations of the mean H I Lya absorption. 
The curves represent He I 584 A absorption for the cases (1) an 
IGM with aiHell = 1 (thick long-dashed curve with pluses), (2) 
the same but masking regions that correspond to ^1,1216 > 3 in 
the coeval Lya forest (thin long-dashed with x's), (3) the He I 
584 A line using our model with THell = 10 -16 s — 1 in the voids 
(thick solid with filled squares), and (4) the same but again masking 
regions that have Thi 1216 > 3 and multiplying by 10 in both panels 
(thin solid with open squares). Also shown are T e fj and a for the 
foreground Lya absorption that falls on top of the 584 A signal 
from 2584 (thick short-dashed curve), and the same but masking 
regions with r > 0.2 to suppress this Lya foreground (thin short- 
dashed curve). 

(|Furlanetto fc Oh! 120081 : IMcQuinn etal] 12009). Even in 
He in regions, the densest gas parcels would have self- 
shielded such that the helium remains largely singly ion- 
ized. Therefore, a detection of He 1 absorption will not 
necessarily indicate that He 11 reionization was occurring, 
but a detection of He 1 absorption from low density gas 
would. The lower the density in which this absorption is 
detected, the stronger the constraint on rueii (and £Heii) 
in models where the gas was ionized by a weak back- 
ground. Furthermore, the H 1 Lyman forest provides 
a measure of the density of an absorber, which allows 
one to target low density gas parcels. We showed in Sec- 
tion that the amount of He 1 absorption in regions with 
7111,1216 < 3 (which roughly correspond to regions with 
Aft < 1) is drastically different between such models. 

If absorption from low density gas were detected, it 
would need to be proven that it is not due to a low-z 
interloper in order to claim a detection of a large-scale 
He 11 region. This section discusses two methods to es- 
tablish that He 11 reionization was occurring at z ~ 3: 
(1) a direct method that involves identifying individual 
weak He 1 lines, and (2) a statistical method that uses 
the cross correlation with the Lya forest to detect the 
He 1 absorption. Unlike the former method, the latter 
method can be performed even when the signal-to-noise 
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ratio (SNR) on individual He I lines is much less than 
unity. 

4.1. Direct Identification 

With a high enough SNR per resolution element, it 
may be possible to dir ectly detect indiv i dual H e I ab- 
sorbers as was done in iReimers &: Vogell (|1993[ ) (there, 
for relatively dense systems). Once a candidate absorber 
is identified, it then needs to be ruled out that the ab- 
sorption is not from a foreground interloper. To detect 
He II gas from regions near the cosmic mean density 
requires targeting He I absorbers with THei,584 ~ 0.1 
(thi,1216 ~ 3). Here is an example of how direct identi- 
fication might be performed: For wavelengths that cor- 
respond to He I absorption at z ~ 3, we count w 10 H I 
Lya systems in 10 4 km s _1 of spectrum with a maximum 
flux decrement between 5 and 10 percent. The probabil- 
ity that one of these systems falls within 50 km s _1 of 
a single He I absorber with the same decrement is 10% 
and of 2 He I absorbers is 1%. Thus, by targeting He I 
absorption from coeval H I systems, one can rule out 
chance contamination from foreground Lya once a few 
systems are detected. 

An alternative method for direct identification would 
be to find a second He I line (such as the 537 A line) 
from the same absorption system. However, because the 
absorption in other He I lines will be at least several 
times weaker, it will be difficult to detect a second line 
from systems near the cosmic mean density (requiring 
SNR ~ 100 per spectral pixel). 

4.2. Cross Correlation 

It is most likely that the SNR per resolution element 
will not be high enough to directly identify He I lines 
in regions with A& ~ 1 with present instruments. Since 
THei,584 °c x He ii thi,1216) with < x H cii < 1 and poten- 
tially more structure in 7^1,1216 than in XHcii, a mea- 
surement of thi,i216 may be used to construct a model 
for the He 1 signal that can be used to extract it in cross 
correlation. In what follows, we denote the He 1 trans- 
mission field as TkcLi, where i labels the spectral pixel. 
Furthermore, we assume Fhoi/Fhi = 1 and XHeii = 1 
to estimate the He 1 transmission from thi,1216- Thus, 
our estimate for the He 1 transmission in a pixel is 
7jfli,i = exp(-0. 025 thi, me,*) ( cf - equations H and 0}. 
However, in the limit TH e l,584 ^ F the SNR of the cross 
correlation is not affected by what we assume for the 
proportionality between THei,584 and ihi,i216- In fact, we 
find that it is only marginally sensitive to this propor- 
tionality in realistic cases because the He 1 absorption is 
weak and because the SNR in cross correlation depends 
more on the phase of modes than on their amplitude. We 
discuss in Sec. 14.2.21 how an imperfect estimate for T£l\ 
affects our estimates. 

An estimator for detecting this signal in cross correla- 
tion is 



turn* fH C i(k)\ v&(k)*D(k) 



\T£im\ 2 



Pm{k)+P ins t' 



(6) 



where D represents the data, tildes denote a quantity 
in Fourier space (a convenient basis because the covari- 
ance matrix of the noise is diagonal), Phi is the power 
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Fig. 4. — Average SNR for detection of He I absorption in cross 
correlation. This quantity is plotted as a function of <7i ns t (the in- 
verse of the average SNR on the continuum in a resolution element), 
assuming a spectrum with TZ = 20, 000 and that the measurement 
uses 1000 resolution elements (Az ~ 0.2). In the top panel, the 
thick curves correspond to z = 4, and the thin to z = 3. The 
solid curves are the full He I signal, and the lower sets of curves 
mask overdense regions that have coeval H I Lyo optical depths of 
thi, 1216 > 3. The bottom panel investigates the effect of different 
masks based on coeval H I absorption or on the observed optical 
depth t in the 584 A forest. Also shown is the SNR of detection at 
1Z = 2000 and for 1000 resolution elements, again masking regions 
with coeval Tfn^ie > 3. 

spectrum of the foreground H I absorption, and Pi ns t is 
the power spectrum of the instrumental noise (which we 
assume is white). Conveniently, the variance on this esti- 
mator owing to noise is its average, making it an estimate 
for the square of the SNR. This estimator is unbiased in 
the sense that components that do not correlate with T£l\ 
do not contribute. This estimator assumes that the sig- 
nal is known since it requires 7h c i as input. In practice, 
this assumption can be avoided because for interesting 
cases we find that T£l\ oc 7h i (and the normalization 
factor can be determined in a measurement to fractional 
precision » 0.5 SNR" 1 ). 

The ensemble average of this estimator over the noise 
and foreground absorption yields 



Vfe 



r(fc) 2 |7Hci(fc)| s 
P H i(fc)+Pi 
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The cross correlation coefficient r(k) is bounded such 
that —1 < r(k) < 1 and is defined as 



r(k) 



75g*(*QT H ciW 

(I^WH^WI 2 ) 1 / 2 ' 



(9) 



The approximate equality given by equation ((5J) is a 
rough estimate for the SNR of detection if r w 1 and 
if the information is coming from modes near the reso- 
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lution limit (which roughly holds for TZ < 20,000). In 
this equality, A^es is the number of resolution elements, 
and 0H eI , cr^j and of nst are respectively the variance 
of the He I 584 A transmission (normalized such that 
T"Hei,584 = is transmission of unity), the foreground 
H I transmission, and the instrumental noise in a reso- 
lution element. This approximate relation shows that, 
for TZ = 2 x 10 4 and N res = 10 3 (Az ks 0.2), such a 
cross correlation can be used to detect the signal even if 
the fluctuations in the signal are more than an order of 
magnitude smaller than the fluctuations in the noiseQ 

Figure @] plots the expected average SNR as a function 
of Cinst at which the He I absorption can be detected in 
cross correlation, assuming N ICS = 10 3 and TZ = 2 x 10 4 . 
To calculate this quantity, we substitute terms of the 
form TT with their ensemble-average in equation ([?])• 
The red curves are this SNR for observations of He I 
absorption at z = 3 (thin) and z = 4 (thick). He I 
absorption can be detected statistically using this tech- 
nique even with an extremely poor SNR per resolution 
element (at » 10 a for a- lns t = 0.3). In the case of COS, 
N IeB « 500, 1000, and 1500, depending on whether 1,2, 
or 3 of the w 40 A bandpasses (that each COS NUV 
grating can simultaneously observe) cover the He I forest 
for a given target. 

For the calculations in Figure EJ wavevectors only up 
to the Nyquist cutoff for the spectrograph resolution are 
included (/cn = ttTZ/c) , with this cutoff imposed using a 
tophat hlter in fc-space. Spectrographs typically sample 
the signal S(n) = J d\' W n {n AA - A') «F(A') to wavevec- 
tors that are a few times higher than k^ to avoid aliasing 
(i.e., A/ A A ~ few TV), where n is the integer number of 
the spectral pixel, T is the incident flux, and Wn is the 
detector window function that damps wavevectors above 
fcN. For realistic a lnst and TZ > 20,000, most contribu- 
tions to the He I SNR of detection are coming from k- 
modes a couple times smaller than the Nyquist wavevec- 
tor (from wavevectors where Pi ns t ~ -Phi), such that our 
results are somewhat robust to the exact details of the 
spectrograph. This can be noted in Figure [5j which fea- 
tures the power spectrum of the different components of 
the signal. The rightmost vertical line is the Nyquist 
cutoff for TZ = 20, 000. The foreground Lya power spec- 
trum (dashed green curve) intersects the noise power for 
0mst ~ 0.05 and <7i nst ~ 0.1 (dotted black diagonal lines) 
prior to the Nyquist cutoff. 

This point of intersection also shows that both the 
instrumental noise and the noise from foreground ab- 
sorption are of comparable importance in determining 
the SNR. Interestingly, because of the different ther- 
mal broadening scales, the He I absorption has signif- 
icantly more power than the foreground absorption at 
k > 0.1 s km -1 (cyan dot-dashed curve, Fig. [S]). This 
difference increases the detectability of the He I absorp- 

7 The discussion in this section could equivalently be phrased 
in terms of optimal signal extraction with a matched filter, where 
the matched filter is Mi = [C^ 1 ]^ 73?j • and Cjv is the covariance 
matrix of the noise (composed of both the instrumental noise and 
low-z Lya absorption). For 7||Ji j = 7hcIJ> Mi is the filter that, 
when its dot product is taken with the signal, maximizes in the 
filtered data the ratio of power in the signal to the average power 
in the noise. Furthermore, equation J7J) is the SNR that such a 
filter can achieve. 



tion in TZ > 20, 000 spectra. The H I foreground absorp- 
tion is also highly non-Gaussian, which could affect the 
probability distribution of the estimated SNR from equa- 
tion ([6]), especially for short skewers that do not sample 
a representative set of foreground absorption lines. We 
have calculated the probability distribution of the esti- 
mated SNR from mock observations constructed from 
our simulations. In practice we find that the probability 
distribution is relatively Gaussian even for 25 /h cMpc 
skewers. 

4.2.1. Masking to measure He I absorption in low density 

regions 

To isolate the He I absorption from low density gas, we 
mask pixels in the simulated He I forest with coeval Lya 
absorption of thi,i216 > 3. This mask covers 10% of the 
pixels at z = 3 and 25% at z = 4. The calculation of the 
SNR in the presence of a mask is more complicated and 
is no longer diagonal in Fourier space. In addition, the 
mask artificially correlates noise wavevectors with those 
of the signal. To avoid such issues, our calculations only 
mask the signal (by setting thi,1216 = 0) an d not the 
noise to calculate the SNR via equation ©. This esti- 
mate for the SNR likely underestimates the value that an 
optimal analysis could achieve because it involves adding 
additional noise to the signal. (See Appendix 151 for ad- 
ditional discussion with regard to the effect of a mask.) 

The full calculation in the presence of masks is com- 
putationally tractable (with just the additional compli- 
cation of inverting Afpixei x Afpixei matrixes, where Npix e \ 
is the number of spectral pixels) . In addition, a simpli- 
fied analysis that applies the cross-correlation statistic 
(eqn. [H]) to just the intervals between masked regions, 
and sums the SNR 2 from each of these intervals, would 
likely provide a reasonable estimate for the SNR of de- 
tection. 

For the mask with thi,i216 > 3, the standard deviation 
of the He I absorption is reduced by approximately a 
factor of 10 for the model with cchcII = 1 relative to the 
unmasked model, and by a factor of more than 100 for the 
model with ThgII = 10~ 16 s _1 (Fig. [3J. However, even 
when masking thi,1216 > 3, it is still possible to detect 
the signal in the model with £Heii = 1 with N ICS ~ 10 3 for 
Cinst = 0.1 (the dashed green curve in the bottom panel 
of Fig. 0]). This value <7i ns t = 0.1 may be achievable with 
COS on several existing sightlines, as discussed in the 
introduction. 

Thus far we have discussed masking regions that corre- 
spond to tjji.1216 > 3. Such a mask requires a SNR w 50 
on the H I Lya forest continuum at TZ > 2 x 10 4 to be 
able to detect flux in pixels with thi,i216 = 3 at 2 a. Such 
a SNR may not be achievable on all He I forest quasars. 
If instead only SNR w 20 can be achieved, then mask- 
ing thi,i216 > 2 is more reasonable. The signal can still 
be detected at 3 a with such a mask with <Ti nst s» 0.06, 
TZ = 20, 000, and AT rcs = 1000, as illustrated by the red 
solid curve in the bottom panel of Figure |4] In addi- 
tion, if Ly/3 information is used, masking only regions 
with thi,i216 > 10 is even possible. Such a mask results 
in much more signal for the £Hcii = 1 model (orange 
dot-dashed curves) , with a 3 a detection possible even at 
Cinst ~ 0.3. Even though this mask includes denser re- 
gions than the mask with thi,1216 > 3, the expected SNR 
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Fig. 5. — Square root of the dimcnsionless absorption power spec- 
trum. The dot-dashed cyan curve is this for the He I absorption at 
2 = 4, and the thick dashed green is this for the foreground Ly« 
absorption. The solid red curve is this for the 2 = 4 He I signal after 
masking regions with coeval thi,1216 > 3 (the different scaling at 
high-fc compared to the blue curve owes to the mask). The vertical 
lines are the Nyquist cutoff for TZ = 2000 and TZ = 20, 000. The 
black dotted diagonal lines are the instrumental noise assuming 
TZ = 20, 000 and either cr inst = 0.05 or a inst = 0.1. 

is a factor of 3 times larger at cr ms t ~ 0.1 in the xudi = 1 
model than in the model with ThcII = 1CP 16 s _1 @ 

4.2.2. Imperfect signal template estimates 

The exact form of the He I absorption signal will never 
be known, and instead we have some best guess T^l\ i for 
its form from coeval Lya forest measurements. In cases 
where the signal is not perfectly known, the achievable 
SNR is decreased by the cross correlation coefficient r 
(eqn. [7]). The different amounts of thermal broadening 
between He I and H I absorption lines make r < 1 even 
with a perfect measurement of thi,i216- This effect was 
included in all our calculations, but we find it has a neg- 
ligible effect on the SNR for detection of the cross cor- 
relation (< 10%). To the extent the IGM had a single 
temperature, the effect of thermal broadening is a convo- 
lution in wavelength space with a single Gaussian filter. 
In this single-temperature limit, the filters divide out and 
do not affect the value of r and, therefore, the SNR of 
detection. 

Another complication arises because the measurement 
of thi.1216 will never be perfect. There will be noise in 
the estimate for this field that reduces r further. To test 
the importance of an imperfect reconstruction of thi,i216j 
we have added noise to thi,i216 with the standard devi- 
ation in the noise per pixel of cr™ t /exp(— Tai.me)- We 
find a negligible difference in the SNR of detection for 
<T hast ~ 0.1 for the case with the thi,i216 > 3 mask, which 

8 The data can be coaxed further to improve the significance 
of detection. For example, much of the contamination from fore- 
ground H I absorption can be eliminated by removing all pixels that 
have t > 0.2 (see the thin green dashed curves in Fig. [3j|. The 
motivation for this additional mask is that the foreground absorp- 
tion is characterized by rare absorbers that typically have larger 
optical depths than the He I absorbers. This operation improves 
the SNR further (blue dotted curve in the bottom panel of Fig. 
2J. Although, such an operation would make the full analysis more 
complicated because it correlates the foregrounds with the mask. 
In addition, masking lines based on their associated Lyman-series 
and other metal lines would further improve the SNR. 



indicates that noise in the reconstruction of Tni.1216 does 
not significantly degrade the measurement 

The He I absorption signal can also differ from the es- 
timate for this signal from thi,i216 because of the patchy 
structure of He n reionization. The signal in the cross 
correlation is coming from k ~ 1 cMpc -1 , whereas the 
structure in the He n fraction is thought to be modu- 
lated at scales of i?bubbio ~ 30 cMpc, the size of the 
He hi bubbles around quasars (|Furlanetto k. Dixonll2009l ; 
iMcQuinn et al.l r2009f ). The wavevectors affected by the 

structure of He n reionization k ~ -^bubble W H1 niore 
than an order of magnitude smaller than those that 
contribute to the SNR. If Hell reionization is perfectly 
patchy such that the helium is cither singly ionized or 
doubly ionized and a skewer of length 3> i?bubbie is mea- 

_ 1/2 1 

sured, then r « x HoII at k ^> "bubble' where Xn e ii is 
the volume-averaged He II fraction. Therefore, the total 
SNR of the cross correlation is reduced on average by the 
factor XHeii- Also note that in the contrasting toy case 
that He n reionization was perfectly homogeneous, the 
SNR is also reduced by XHcii- 

A detection of He I absorption from low density gas 
would translate to a Iff error on ShoII of 5xn c u ~ 
xfj eII (SNR)/SNR, where x| cII (SNR) is the expected 
a^Heii for a given estimated SNR of detection in cross 
correlation with exp(— 0.025 thi,i216)- Thus, a 4cr detec- 
tion from low density gas would yield a 25% constraint 
(ignoring uncertainty in Thci)- This constraint assumes 
that the low density gas during He n reionization was 
composed of large-scale regions in which the helium was 
either He n or He III, which is expected theoretically and 
is what is seen in simulations (McQuinn ct al. 2009). 

5. CONCLUSIONS 

This paper discussed the usefulness of the He I 584 A 
forest to study He II reionization at 3 < z < 4.5. 
The optical depth of the He I 584 A line is propor- 
tional to the optical depth in H I Lya by the factor 
0.025 XHoii Thi/ThcI (ignoring differences in the amount 
of thermal broadening between H I and He i). The factor 
XHoii should have been spatially variable, but the ratio 
Thi/ThcI should have been essentially spatially indepen- 
dent and roughly equal to unity at relevant redshifts. 
Therefore, this absorption can be used to study He II 
reionization through its dependence on XHeii- 

Our best method at present to probe He n reionization, 
He II Lya absorption, saturates at neutral fractions of a 
part in a thousand at the cosmic mean density. In con- 
trast, He I 584 A absorption is unsaturated in all except 
the densest regions, even for XHcii = 1- This absorption 
provides a complementary window into the ionization 

9 We found earlier that noise at the 10%-level in the NUV spec- 
trum can have a substantial effect on the SNR at which the He I 
absorption is detected. The reason that a similar level of noise is 
less important in the H I spectrum is because tjji 1216 ~ 1 f° r most 
pixels at these redshifts, which translates into a ~ 10% error on 
T HI,1216 if °hist = except in the low optical depth pixels (where 
there is little He I signal anyway) and the high ones. Whereas for 
He I absorption, if <?P|J: = 0.1, we have order unity uncertainty in 
this optical depth for typical THel,584- Furthermore, r (and thus 
the SNR of detection) is most sensitive to the phase of modes, 
which is largely maintained even when significant noise is added to 
the high- 7hi, 1216 pixels. 
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state of intergalactic helium at z ~ 3 that can defini- 
tively test whether an opaque region in the He n Lya 
forest was due to a large-scale He n region. Even for re- 
alistic amounts of instrumental noise and foreground H I 
Lya absorption, we showed that the coeval H I Lya for- 
est absorption can be used to construct a matched filter 
that can detect the He I absorption at high significance 
from a single quasar sightline. 

A detection of He I absorption from A& ~ 1 gas at 
z ~ 3 would definitively indicate that He II rcionization 
was occurring. If a significant fraction of the intergalac- 
tic helium was in He II, we found that He I 584 A forest 
absorption from low density gas could be identified in 
a quasar spectrum with SNR ~ 10 and TZ ~ 10 4 in 



an interval of Az = 0.2. These specifications may be 
achievable with the COS instrument on the HST for a 
few known z > 3 targets. 
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APPENDIX 

A. MEASURING THE HARDNESS OF THE IONIZING BACKGROUND WITH He i ABSORPTION 



We argued in Section lPl that the spectrum between 13.6 eV and He II Lya at 40.8 eV (and, thus, the ratio Thi/Thci) 
depends only on the H I column density distribution and the spectrum of the sources. Since the H I column density 
distribution i s constrained by observations, a measurement of Fht/ ThrT will constrain the spectrum of the sources 
in this band (jMiralda-Escude fc Ostrike"r1ll992t ISantos h Loebl [2003 ) . This ratio can be measured by comparing the 
conti nuum absorption at 912 (1 + z) A - the H I limit - and 504 (1 + z) A - the He I limit (Miralda -Escude fc Ostriker] 
1992). These two optical depths are related via 

cont l pion \ _ „ ""He a HcI I'm 'Hcl(-5'Hcl) ront I P HoI ^ 

-Hcl [E m ) ~ XHel! — — ^ ^ion) %I ). 

wO-lOXHell^TST ( Al ) 
1 Hel 

where primes denote the average value for the quantity inside the absorber. This ratio of continuum optical depths 
can be measured most easily in systems with JVhi ~ 10 17 cm~ 2 , corresponding to systems with r^ nt (E^) ~ 1. These 
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column densities will be extremely self-shielded to He n-ionizing photons such that XHeii = 1 is a good assumption. 
Thus, the ratio of the continuum optical depths can be used to measure rj iI /r^ icI . Because this method relies on 
detecting smooth spectral features, this science can be pursued even in much lower resolution spectra than are required 
to study He n reionization, and a < 10% constraint on the continuum is required to detect the weak breaks from He I 
continuum absorption in systems with T^ nt (E^) ~ 1 and, thereby, constrain r^j/T^j. 

The primary complication with this continuum-absorption method to constrain Tni/THei is that He I self-shields 
at H I column densities that are a decade larger than those for which the H I self-shields. Thus, in systems with 
Th^C-EhS) ~ 1) rjjj < Thi even though r^ cI rs ThcI- However, if the absorber is a slab of material for which both the 
line of sight and the incident radiation is perpendicular to the plane of the slab, then T' m = Thi [1 — exp(— T^° nt )]/T^j nt . 
The difference in this case between and Thi is relatively small for rg° nt < 1. (In fact, under most circumstances, the 
correction will be smaller than in this toy example because there will be other orientations with smaller optical depths to 
locations within the absorber.) Therefore, a measurement of ^H°rV r Hi nt (which yields T^j/r^j) for TVhi ~ 10 cm~ 2 
can be mapped to a constraint on Tni/rud without significant error. 

A second method to constrain Tni/rud is to measure the equivalent width of a l ine in the He I forest a nd to compare 
this with the equivalent width of this system's H I Lyman-series absorption lines (|Santos fc Lo cb 2003). This method 
can be employed on systems that have much smaller H I column densities than the above continuum absorption 
method. What this method requires is absorption systems that have xhoII ~ 1 in order to directly measure Tni/ruei 
(equation|31). For TroII = 10~ 15 s" 1 , systems w ith JVhi > 10 15 cm~ 2 at relevant z self-shield to He Ii-ionizing photons 
and can have ccHeii ~ 1 (jMcQuinn et al.|[2009() . At lower Fueii, even smaller H I column densities self shield. (Note 
that, for a fixed ThcII, the He II is likely to transition from ionized to neutral over a relatively narrow range in TVhi- 
The same does not happen to the hydrogen once it starts to self-shield, and it is like to remain highly ionized for 
Nm < 10 19 cm" 2 at z ~ 3.) 

The major difficulty with the second method is that it will be challenging to measure the H I column density 
for systems that have a significant optical depth in He I 584 A. The use of higher H I Lyman series resonances, or 
continuum absorption for the largest column densities, will be required for such measurements. For example, systems 
with ./Vhi ~ 10 15 cm -2 and XHcii = 1 have THei,584 ~ 1 an d H I Ly7 absorption of thi,973 ~ 2, assuming a linewidth of 
30 km s _1 . Even if XHeii = 1 is not assumed, such a measurement would place an upper bound on Thi/FucI- 

B. SIGNAL-TO-NOISE OF CROSS CORRELATION 

Denote the measurement and estimated templates with subscript m and e, respectively. Both the template and 
measurement have some underlying signal and noise, which we represent here for the measurement as 7~ m {k) — 
T m ,t{k)+n m {k), where T m ,t is the underlying He I signal and n m (k) is the measurement noise. A similar decomposition 
applies for the estimated signal: T e (k) — T e ,t(k) +n e (k), where T e ,t(k) is the component of the estimate that correlates 
with the signal and n e (k) is the component that does not. The covariance of the cross-power between the He I 
measurement and the template estimate, H m ^ e (k), is then 

{5H rn Aki)5H* m ^k 2 )) = (H m , e (ki)H* m Jk 2 )) - (H m , e (ki))(H^ e (k 2 )), 

= (T m (k 1 )T;(k 1 )TZ(k 2 )T e (k 2 )) - (T m (fc 1 )T e *(fci))(r^(fc 2 )r e (fc 2 )) ] 

= <[T m , t (fci) +n m (k 1 )][T* t (k 1 ) +n* e (k 1 )][T^ t (k 2 ) +<„(fc 2 )][r e , i (fc 2 ) + n e (fc 2 )]>, 

-{[T m , t (h) + M*i)][7£(fci) + <(fci)])([r^ t (fc 2 ) + n* m (k2)][T e>t (k 2 ) + n e {k 2 )]), 
= 7^,t(fci)7^* !t (fc 2 )(n*(fci)n e (fc 2 )) + 7^* t (fci)7^,t(fc 2 )(n m (fci)n^(fc 2 )) + (n m (fci)n^(fc 2 ))(n*(fci)n e (fc 2 )), 
= 4 1 ,fe 2 [|r m , t (fci)| 2 A^ e (fci) + ITeAh^N^h) + N m (ki)N e (ki)], (Bl) 

where Nx(k) = (nx{k)n* x {k)) and in general, Tx,t{ki)T x t {k 2 ) is non-diagonal, but is multiplied here by 
{rL* x {ki)nx{k 2 )) 1 which is 8k lt k 2 ^x{ki). The significance with which we can measure the amplitude of the cross- 
power is 

snr 2 = J2 r e , t (fci)r,:^(fc 1 )( ( 5i? m ^(fc 1 ) ( 5 J ff m , e (fc 2 ))- 1 r e : 4 (fc 2 )r„ M (fc 2 ), 



k± ,k,2 



E 



\TeAk)T* t (k)\< 



. [\TrnAk)\ 2 Ne{k) + \Te,t (*1 ) 1 2 N m (fc) + N m (k)N e (k)] ' 

v \%Ak)T^ t (k)\ 2 

~2-,\%Aki)\ 2 N m {ky [ ' 

where the last line assumed that template estimate noise N e (k) is negligible. Ignoring the additional He I ionization 
structure from fluctuations in XHeii (which we argued in Section [4.21 mostly affects the amplitude of small-scale modes 
in Tm,t), one has N m (k) = -Pmst(fc) + Pm(k). With this replacement, Eg. IB2I coincides with Eq. [71 

In the text, we estimated the SNR in the case where only the signal is masked but not the noise. The derivation 
of Eq. IB2I holds in this case, except that 7^,t(fc)7^ t (fc) is replaced with this cross correlation between the (noiseless) 
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estimated and measured signal, including the masked regions. Furthermore, one could imagine adding mock noise in 
masked regions to real observations to replicate this case. 

In a proper analysis of observations, the noise will be masked along with the signal. Since the mask covers a relatively 
small fraction of space at relevant redshifts, the estimate for the SNR given by Eq. lB2l should approximately hold. This 
can be noted by first considering the SN R on e can achieve if the segments between masks were considered separately. 
In this case, the SNR 2 is the sum of Eq. IB2I computed from each segment. The difference in this case is (1) that the 
SNR is decreased by « (1— /mask) 1 / 2 where / mas k is the mask covering fraction because there are / mas k fewer k values, 
and (2) T e ,t(k)Tm t (k) in Eq. IB2I is replaced with its value in the unmasked regions (which is slightly different than the 
replacement discussed in the previous paragraph, which included the power from the mask). This argument ignores 
the correlations that occur between the noise from the foreground H I absorption in different segments. However, these 
correlations should be small since the low-z forest is characterized by discrete, largely uncorrelated lines. 



